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Abstract 

During the past years several variance reduction techniques for Monte Carlo electron transport have been developed 
in order to reduce the electron computation time transport for absorbed dose distribution. We have implemented 
the Macro Response Monte Carlo (MRMC) method to evaluate the electron spectrum which can be used used as a 
phase space input for others simulation programs. Such technique uses probability distributions for electron histories 
previously simulated in spheres (called kugels). These probabilities are used to sample the primary electron final state, 
as well as the creation secondary electrons and photons. We have compared the MRMC electron spectra simulated in 
homogeneous phantom against the Geant4 spectra. The results showed an agreement better than 6% in the spectra 
peak energies and that MRMC code is up to 12 time faster than Geant4 simulations. 
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1. Introduction 

Currently, modeling electron beam treatment 
planning system is a challenging task. Several 
methods have been provided to calculate the tridi- 
mentional (3D) dose distributions, from Hogstron 
analytical approximation [1] to Monte Carlo (MC) 
simulations. The simulations with MC code, using 
condensed history, are accurate but usually require 
several hours to achieve reasonable statistics. To re- 
duce the simulation time several variance reduction 
techniques have been proposed such as Response 
History Monte Carlo (RHMC) [2j, Macro Monte 
Carlo (MMC) 01], Phase-space evolution Monte 
Carlo (Evolution) 0, Macro Response Monte Carlo 
(MRMC) and Voxel Monte Carlo (VMC) 0. 
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The RHMC, MMC, Evolution and MRMC codes 
use probability distributions generated from elec- 
tron histories previously simulated in a volume ele- 
ment, to transport the primary electron and create 
secondary ones. The Evolution generates its phase 
space in cubic voxels, the RHMC generates in hemi- 
spheres and MMC and MRMC generate in spheres 
(kugels). The VMC code relies on some simplifica- 
tions of the models for electron transportation and 
history repetition. 

The cited codes have shown an accuracy up to 
5% in the depth dose, but only RHMC, MMC and 
VMC have achieved a real speed gain, when com- 
pared with others MC codes. Basically, the use of 
voxels instead of kugels is the limitation of the Evo- 
lution code because it is difficult to overcome the 
symmetry problem. The MRMC limitation, men- 
tioned by the author, was related to the small radii 
of the kugels used. 

We have implemented the MRMC method with 
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larger radii as well more radii choice for electron 
transportation. Such implementation can be used 
for generating electron spectra, since it does not have 
dose deposition algorithm implemented. 

2. Methods 

2.1. MRMC implementation 

The MRMC simulation is divided into two parts 
called local and global calculations. In the first one 
we have simulated the kugels database for water, 
soft tissue and compact bone. In the second part is 
the electron interaction simulation in the materials 
and phantom geometries selected. 

2.1.1. Local calculation 

The local calculation consisted in the simulation 
of electron histories starting from kugel's center with 
initial direction aligned to the Z axis, figure [TJ The 
particles' state leaving the kugel is tailed in proba- 
bility distribution histograms. 




Figure 1. Illustration of kugel bands with some electron 
tracks and its output band. 

The kugel is divided into 4 bands along the Z axis, 
where the projections over such axis are: 1.0 to 0.98 
(band 1), 0.98 to 0.8 (band 2), 0.8 to 0.4 (band 3) 
and 0.4 to -1.0 (band 4). Each of these bands has a 
complete set of probabilities distribution histograms 
comprising one for the output positions, two for di- 
rections and one for energies. These quantities have 
been divided into three groups: primary electrons, 
secondary electrons and photons. Also, a histogram 
had been generated to store the average number of 



escaping secondary electrons from the kugel geom- 
etry and another one for photons escaping. 

Each material has a complete set of several kugel 
radii and energies. The radii chosen were: 0.025, 
0.050, 0.100, 0.250, 0.500 and 1.000 cm, where only 
the first three were tested in the original work. We 
have chose 19 energies stepping from 31 MeV down 
to 178 keV with an energy spacing around 25%, as 
suggested by Svatos [6] . 

We have used the Geant4 version 8.0 [8.9J to simu- 
late the kugels database. Such code uses a condensed 
history for electron transportation, while in the orig- 
inal implementation of the MRMC the author have 
used a Single Scattering MC code. Each kugel has 
been simulated with 5 x 10 5 histories, using 990 eV 
as electron cutoff and 10 ^m as rangecut for water 
and soft tissue and 20 ^tm for compact bone. Due to 
long elapsed time, these simulations took place in a 
small cluster. 

The database was stored in ROOT [10] file for- 
mat, which properly store kugel data for global cal- 
culations. 

2.1.2. Global calculation 

The global calculation is the electron transport 
through the phantom geometry using kugels to sam- 
ple the primary electron final state, as well the sec- 
ondary electrons and photons production. An activ- 
ity diagram, in Unified Modeling Language (UML), 
is presented in the figure[2j First, the system tries to 
find an energy through a linear-logarithmic statis- 
tical interpolation, expression [H where E is the en- 
ergy of electron being transported, E+ and E- are 
the nearest energies available in the database and x 
is an uniform random number in [0, 1[. The E + is 
chosen when the expression Q] is true, otherwise E_ 
will be chosen. 

logE+ - logE 
X logE+ - logE- [ ' 

After an energy is chosen, the system searches for 
a radius that does not cross any boundary. If there 
is a kugel, the MRMC code system uses it to sample 
the particles final state, including the secondaries 
and photons production. If there is no kugel, the sys- 
tem transports the current electron using an straight 
line approximation and corrects its energy using a 
Continuous Slowing Down Approximation (CSDA). 
The secondaries are transported as the primary ones 
and the photons are transported until they leave 
the geometry, since interaction models are no imple- 
mented. 
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Figure 2. UML Activity diagram showing the electron logic 
transportation. 

The MRMC was written in object oriented style 
using C++ language, with the ROOT library to 
sample the histograms and for the particle trans- 
portation in the geometry. 



agreement better than 2% in the peak energies for 
the incident energies, above 3 MeV, and up to 5% 
for incident energies below 3 MeV, figures [3] to [U 
For compact bone the results showed an agreement 
better than 6% for the incident energies above 
5 MeV and better than 2% for the incident ener- 
gies below 5 MeV, figures [6] and [Jj Such deviations 
may be due to uncertainties in the local simulations 
done with Geant4, straight line approximation for 
boundary cross and uncertainties in total stopping 
power tables [H] . 
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Figure 3. Forward electron spectrum for incident energy of 
10 MeV in water phantom. 



2.2. Benchmarks 

The tests consisted in an electron pencil beam in- 
cident perpendicularly in a homogeneous cylindrical 
phantom. The electron beam starts at 0.025 cm, in- 
side the phantom, in order to avoid the straight line 
approximation in the beginning of event. The simu- 
lation parameters for the MRMC were set to 5 keV 
for both, the cutoff and secondaries threshold. The 
Geant4 parameters were set for 1 keV of cutoff and 
0.1 mm of rangecut. 

We compared the forward and lateral electron 
spectra scatter against the Geant4 simulations. The 
energies chosen for these tests were: 10.0, 8.0, 7.5, 

5.0, 3.0 and 1.0 MeV. 

3. Results 

3.1. Electron spectra 

The forward and lateral spectra simulated by 
MRMC for water and soft tissue have shown an 




Figure 4. Lateral electron spectrum for incident energy of 
10 MeV in water phantom. 



3.2. Simulation speed 

The MRMC simulation speed was around three 
times faster for lowest energy and up to 12 times 
faster for 10 MeV. Simulations for low energies are 
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4. Conclusions 
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The simulated spectra from our implementation 
of MRMC code have shown good agreement and the 
speed was up to 12 time faster, when compared to 
Geant4 for energies from 1 to 10 MeV. 

The authors would like to thank to prof. Dr. An- 
tonio Carlos Roque da Silva Filho and Dr. Rodrigo 
Freire Oliveira to allow the access to cluster for our 
calculations. Also, we would like to thank to CNPq 
for the financial support. 



Figure 5. Forward electron spectrum for incident energy of 
3 MeV in soft tissue phantom. 




Figure 6. Forward electron spectrum for incident energy of 
5 MeV in compact bone phantom. 



Lateral electron spectrum - 5 MeV 
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Figure 7. Lateral electron spectrum for incident energy of 
5 MeV in compact bone phantom. 

slower because the database for such energies have 
only smaller kugels radii available. 



4 



Espectro frontal dos eletrons - 5 MeV 
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Espectro lateral dos eletrons - 5 MeV 
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